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We present a theory to investigate electro-kinetic behavior, namely, electrorotation and dielectrophoresis un- 
der alternating current (AC) applied fields for a pair of touching inhomogeneous colloidal particles and biologi- 
cal cells. These inhomogeneous particles are treated as graded ones with physically motivated model dielectric 
and conductivity profiles. The mutual polarization interaction between the particles yields a change in their 
respective dipole moments, and hence in the AC electrokinetic spectra. The multipolar interactions between 
polarized particles are accurately captured by the multiple images method. In the point-dipole limit, our theory 
reproduces the known results. We find that the multipolar interactions as well as the spatial fluctuations inside 
the particles can affect the AC electrokinetic spectra significantly. 
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I. INTRODUCTION 

Identification and analysis of cell populations and (mi- 
cro )biological particles are essential in many practical ap- 
plications ranging from cancer research to chemical anal- 
ysis of environmental pollutants. During the past decade, 
alternating current (AC) electrokinetic phenomena, and in 
particular electrorotation (ER) and dielectrophoresis (DEP), 
have received much attention in this respect, especially in 
micromanipulation and separation of submicron size parti- 
cles CI 011 BH00ll0[IolGjl0- In addition to 
biological and environmental applications, AC electrokinetic 
phenomena have been suggested as possible mechanisms for 
nanomotors PH . 

Both dielectrophoresis and electrorotation are based on di- 
electric properties of particles. These properties depend heav- 
ily on the nature of the surface, e.g., size, shape, and charge 
density. For example, since the composition and shape of can- 
cer cells differ from those of healthy cells, these difference are 
reflected in their characteristic dielectric properties which can 
be exploited in identifying them. From a practical point of 
view, AC electrokinetic methods have the advantages of short 
detection times and high sensitivity |5J. 

Dielectrophoresis can be defined as the movement of 
polarizable particles in a non-uniform applied AC electric 
field 1 15], whereas in electrorotation an interaction between 
a rotating AC electric field [6] and suspended particles leads 
to a rotational motion of the particles. The most commonly 



used models to deal with the dielectric properties of colloidal 
particles or biological cells are the so called shell models. 
Because of inhomogeneous compartmentalization of biolog- 
ical cells, one-, two-, and three- shell models have been ap- 
plied to discuss electrorotation of biological cells, e.g., see 
Refs. BE0. 

These cell models have several limitations and they become 
complex as the number of shells increases. This is particu- 
larly true when two (or more) particles approach each other. 
In the dilute limit, one can focus on the electrokinetic spec- 
tra of an individual particle. If the suspension is not dilute, 
as it is often the case in practice, the situation is complicated 
by the existence of multipolar interactions. Even when a sus- 
pension is initially in the dilute limit, particles often aggregate 
due to the presence of an external electric field. In this case, a 
point dipole approximation 118l 1 1 911 becomes inadequate and 
the mutual interactions must be taken into account by a the- 
ory l20ll2lll that goes beyond the point dipole. 

To provide a physically motivated and tractable model for 
inhomogeneous particles, such as cells, we have recently stud- 
ied particles with spatial gradients in their structures by intro- 
ducing profiles for the conductivities and dielectric constants 
of the particles, and used differential approximation for the di- 
electric factor i22l l23ll24ll . Here, we extend this work to take 
into account polarization interactions when two particles ap- 
proach each other, and treat both DEP and ER using the same 
theoretical framework. We consider a pair of touching graded 
particles in suspensions. As a result, the mutual polarization 
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The DEDA equation for a spherical graded particle has the 
form 1E2H23I1 
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FIG. 1: ER spectra for an isolated particle and two touching par- 
ticles with R/2a — 2.00, 1.10, 1.03, respectively. Parameters: 
c = — 30eo, m — 1.0. 



interactions and the gradation fluctuations inside the particles 
lead to significant changes in the electro-kinetic spectra. In a 
more general context, these are manifestations of correlation 
effects in charge-carrying system I25l l26ll . 



II. FORMALISM 
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It is worth noting that the DEDA is essentially exact since it is 
in an excellent agreement with the exact solutions obtained for 
a power-law profile and a linear profile by solving the Laplace 
equation for the local electric field [ 24] . 

For a pair of particles at a separation R in a suspension, we 
have to consider the multiple images effect |20, 27]. Theoret- 
ically, we may see the inhomogeneous graded particle as an 
effectively homogeneous one. Then, we consider two parti- 
cles in a suspension which is subject to an external uniform 
electric field. This yields a dipole moment into each particle. 
Let us denote the dipole moments of particles 1 and 2 as pio 
and p2o(= P10 for identical particles), respectively. Then, we 
take into account the image effects. The dipole moment pio 
induces an image dipole pu into particle 2, while pu induces 
another image dipole in particle 1 . As a result, multiple im- 
ages are formed. The same description holds for p20- Thus, 
we admit the infinite series of image dipoles. To this end, we 
obtain the sum of dipole moments inside each particle, and 
derive the desired expressions for dipole factors. Let us con- 
sider two basic cases: 1) longitudinal field (L), where the field 
is parallel to the line joining the centers of the particles, and 
2) transverse field (T), where the field is perpendicular to the 
line joining the centers of the particles. 

Based upon a multiple images method Izl . the dipole fac- 
tors, b L * and b T * , are given by fzoHzj 



We consider inhomogeneous biological cells or colloidal 
particles with radius a. We assume that they have a distance- 
dependent complex dielectric constant £i(r) (0 < r < a), and 
that they are embedded in a host fluid having dielectric con- 
stant 62. Here e = e(r) + a(r) / (2mf), where e denotes the 
real dielectric constant, a(r) the conductivity, / the frequency 
of an external field, and i = y— T. As the above formulas sug- 
gest, e(r) and a(r) are not constant inside the particle but have 
distance dependent profiles. This is a very physical assump- 
tion and we will return to it in the following discussion and 
later in connection with the numerical simulations. 

The dipole, or Clausius-Mossotti, factor reflects the polar- 
ization of a particle in a surrounding medium. In a recent 
work 1 22], we derived the dipole factor for graded spherical 
particles by introducing a differential effective dipole approx- 
imation (DEDA). The generalization to the nonspherical case 
was done as well ll23ll . The idea of the DEDA can be sum- 
marized as follows: Consider a shell model for an inhomo- 
geneous particle. In the DEDA one adds new shells of in- 
finitesimal thickness to the particle. Each of these cells have 
distance-dependent complex dielectric constant. Since the 
thickness of the layer approaches zero (d?' — > 0), its correction 
to the dipole factor is infinitesimal and one could eventually 
obtain a differential equation. 
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for longitudinal and transverse field cases, respectively, where 
a satisfies the relation cosh a = R/2a. Although it is not ob- 
vious, it is important to notice that multipoles are included in 
the above formulas l28il . Clearly, the multiple images effects 
have been taken into account in b* L and fry. It is worth noting 
that: Setting n up to 1 in the two equations leads to the dipole 
factor for two touching particles in the point-dipole limit. In 
this case, in view of both \b\ 2 <C 1 and R/2a ~ 1, we have 



™ = TTbTs- (4) 

Both Eqs.@ agree well with the result of Jones, which were 
obtained by a field method in the point-dipole limit 11911 . 
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FIG. 2: ER spectra for two touching particles. Upper panel for profile 
constants m at c = and c = — 30eo, respectively. Lower panel for 
various c at m = and m = 1.0, respectively. The spectrum is given 
as the imaginary part of the dipole factor. Parameter: R/2a=1.03. 



FIG. 3: DEP spectra for two touching particles in the longitudinal 
field case. Upper panel for profile constants m at c = and c = 
— 30eo, respectively. Lower panel for various c at m = and m = 
1.0, respectively. The spectrum is given as the real part of the dipole 
factor. Parameter: R/2a=1.03. 



A. Electrorotation 

By adding a rotating electric field with magnitude E'er* to 
the two particle system, the effective dipole factor for a pair 
of particles should be given by 1 20] 

b* = {b L * +b T *)/2. (5) 

Thus, in this case, the electrorotation velocity of a particle fT 
is given by 1 20] 

Q* = -<Ke 2 ,772,£ER*)Im[6*], (6) 

where 4>( e 2, V2, -E'er) is a function of e-2, the viscosity of 
the medium fy 2 , and -Eer*- Here Im[- • ■] denotes the 
imaginary part of [• •■]. For an isolated spherical particle, 

0(£ 2 ,??2,Eer*) = e 2 E ER * 2 /2r, 2 0. 



B. Dielectrophoresis 

We consider a single particle suspended in a medium and 
subjected to a nonuniform AC electric field E rwp *. The DEP 
force Fdep acting on the particle is given by 1 29] 

F DEP - 27re 2 a 3 Re[6]V|E DEP *| 2 , (7) 

where Edep* stands for the local RMS electric field, and 
Re[- • •] denotes the real part of [■■ •]. Next, for a pair of touch- 
ing particles, the DEP force is given by I2lll2^1 

F L * = 2^e 2a 3 Re[fe L *]V|E DE P*| 2 , (8) 
F T * = 2^e 2a 3 Re[6 T *]V|E D EP*| 2 , (9) 

for longitudinal and transverse field cases, respectively. The 
above formulation for the DEP force is, stricly speaking, ap- 
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FIG. 4: Same as in Fig.|5| but in the transverse field case. 

plicable for linearly polarized fields l30ll . which is the case 
studied here. 



III. NUMERICAL RESULTS 

For the following numerical calculations, we take the con- 
ductivity and dielectric profiles to be 

o-i (r) = <n(0)(r/a) m , r<a (10) 
ei(r) = ci(0) +c(r/a), r<a (11) 

where m and c are profile dependent constants. The profile 
is clearly physical since conductivity can change rapidly near 
the boundary of cell and a power-law profile prevails 1 231 ■ On 
the other hand, the dielectric constant may vary only slightly 
and thus a linear profile suffices |23]. In particular, the dielec- 
tric constant at the center, namely ei(0), may be larger than 
that at the boundary. Thus, in what follows, we would choose 
c < 0. By integrating the dielectric profile, we obtain an aver- 
age dielectric constant e av for different values of c by using a 



FIG. 5: DEP spectra in the longitudinal field case (upper panel) and 
the transverse (lower panel), for an isolated particle, and two touch- 
ing particles with separation ratio R/2a = 2.00, 1.10, 1.03, respec- 
tively. Parameters: c = — 30eo,m = 1.0. 

volume average 1 23] 

£i (r)r 2 dr 

, av _ JO 1 \ I m\ 

ei — r a 9 ; • (12) 

Jo r dr 

For the above dielectric profile t\ av — ei(0) + 3c/4. 

For all numerical calculations, we take ei(0) = 75eo, £2 = 
80e , cri(0) = 2.8 x l(T 2 S/m, and er 2 = 2.8 x 10" 4 S/m. 
Here eo denotes the dielectric constant of the free space. 

Figure shows the ER spectrum of two particles at dif- 
ferent distances from each other. At large separations (e.g. 
R/2a > 2), the multipolar interaction may be neglected, but 
the induced multiple images play an important role in the 
spectrum when two particles approach each other. 

In Fig. [2] it is evident that a second peak due to the mul- 
tiple image effect occurs at a lower frequency. In fact, the 
appearance of a second peak has been predicted for homoge- 
neous particles in a recent work 12011 . Moreover, fluctuations 
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in the conductivity profile can make the characteristic fre- 
quency shifted to lower frequencies (red-shifted), while those 
in the dielectric profile can enhance the peak value. How- 
ever, such effects on the second characteristic frequency and 
its peak value are small enough to be neglected. 

Fluctuations in conductivity and dielectric profiles may en- 
hance the DEP spectrum not only in the longitudinal field 
case (see Fig|3j, but also in the transverse (see Fig@). The 
effects of multiple images may change the DEP spectrum sig- 
nificantly. 

Similar to Fig. [J Figure [5] shows that the multiple images 
play a crucial role in the DEP spectrum when the particles are 
close enough. In contrast to the result from an isolated parti- 
cle, the multiple image effect may enhance the DEP spectra 
(moreover, the real part of the dipole factor may be enhanced 
to be larger than 1) at low frequency region in the longitudi- 
nal field case. However, the DEP spectrum is reduced in the 
transverse field case due to the presence of multiple images. 

In addition, we have also compared the point dipole model 
with the current multiple image dipole model (no figures 
shown). As expected, the results predicted by them are quite 
different, especially at low frequencies. This shows further 
that the point dipole model is inadequate for the touching par- 
ticles, and thus needs to be modified to take into account the 
effect of multiple images. 

IV. DISCUSSION AND CONCLUSION 

In this work, we investigated the effects of multipolar in- 
teractions on AC electrokinetic behavior, electrorotation and 
dielectrophoresis, of inhomogeneous biological cells and col- 
loidal particles. We model such inhomogeneous particles as 
graded ones. Our method may be extended to high concen- 
tration case l3ll l32ll or pearl chain case l33l l34ll . work is in 
progress to address these issues in detail. Also, it is possible 
to take into account shape effects by considering the nonspher- 
ical shapes, such as oblate or prolate spheroid 1 11]. In doing 
so, we might resort to the derived DEDA equation for graded 
spheroidal particles as well J23ll . It is also straightforward to 
extend this work to deal with the experimentally interesting 



case of charged colloidal suspensions 1I2II1 . 

To put our approach in the context of composite particles, 
we have performed a mean-field approach in the spirit of Choy 
et al. 1351 13611 . i.e., treating inhomogeneous particles as effec- 
tively homogeneous ones which are embedded in a uniform 
field. In particular, it is worth noting that well-known Tar- 
tar formula 13711 can be used to exactly calculate the effec- 
tive complex dielectric constant of a single graded particle. 
Thus, once this effective complex dielectric constant is ob- 
tained, one can proceed to calculate the relevant dipole mo- 
ment, and hence the desired dipole factor. More interestingly, 
our DEDA approach can predict exactly the same result as 
Tartar formula Oal . 

To sum up, based on the DEDA, we have presented a the- 
oretical study of electrokinetic behavior, electrorotation (ER) 
and dielectrophoresis (DEP) for two touching inhomogeneous 
particles in suspensions. We found mutual polarization effects 
and the spatial fluctuations inside colloidal particles or bio- 
logical cells can both affect ER and DEP spectra significantly. 
Our approach has the further advantage of being able to treat 
both electrorotation and dielectrophoresis using the same the- 
oretical framework. 

As a further study and a test to our theory, it would be in- 
teresting to have a systematic experimental investigations of 
these effects. The hope is that they would shed light to the 
limits of the theory and that they would help to separate the 
DEP and ER behavior from, e.g. electrohydrodynamic flow 
effects 1 39] and limitations due to Brownian motion. One 
possibility for doing so would be to use the laser tweezers 
combined with ER and/or DEP 
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